Submitted to the Astrophysical Journal 



The Velocity Distribution of the Nearest Interstellar Gas 

Priscilla C. Frisch 
Lauren Grodnicki 
Daniel E. Welty 

University of Chicago, Department of Astronomy and Astrophysics, 5460 S. Ellis Avenue, 

Chicago, IL 60637 

ABSTRACT 

The bulk flow velocity for the cluster of interstellar cloudlets within ~30 pc 
of the Sun is determined from optical and ultraviolet absorption line data, af- 
ter omitting from the sample stars with circumstellar disks or variable emission 
lines and the active variable HR 1099. Ninety-six velocity components towards 
the remaining 60 stars yield a streaming velocity through the local standard of 
rest of -17.0±4.6 km s -1 , with an upstream direction of / n =2.3°, 6 n =-5.2° (us- 
ing Hipparcos values for the solar apex motion). The velocity dispersion of the 
interstellar matter (ISM) within 30 pc is consistent with that of nearby diffuse 
clouds, but present statistics are inadequate to distinguish between a Gaussian 
or exponential distribution about the bulk flow velocity. The upstream direction 
of the bulk flow vector suggests an origin associated with the Loop I supernova 
remnant. Groupings of component velocities by region are seen, indicating re- 
gional departures from the bulk flow velocity or possibly separate clouds. The 
absorption components from the cloudlet feeding ISM into the solar system form 
one of the regional features. The nominal gradient between the velocities of up- 
stream and downstream gas may be an artifact of the Sun's location near the 
edge of the local cloud complex. The Sun may emerge from the surrounding 
gas-patch within several thousand years. 



1. Introduction 



A number of studies have searched for correlations between the velocity of interstellar 
gas observed inside the solar system and towards external stars (Adams & Frisch 1977; 
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McClintock et al. 1978; Lallement & Bertin 1992). Adams & Frisch (1977) showed that the 
velocity of interstellar gas inside the solar system differs by several km s -1 from interstellar 
cloud velocities towards stars located in the upstream direction. The discovery that the H° 
La line is redshifted by several km s" 1 from the D° line towards the nearest star (a Cen, 
1.3 pc) confused the identification of the cloud velocity in this direction (Landsman et al. 
1984). This shift has since been successfully modeled by including the La absorption from 
compressed H° in the heliosheath (Linsky & Wood 1996; Gayley et al. 1997). However, the 
heliosheath is observed only in the strong La line. The D°, Mg + , and Fe + line velocities 
towards a Cen (~-19 km s -1 , Lallement et al. 1995; Linsky & Wood 1996) disagree by ~2 
km s _1 with interstellar gas and dust velocities found inside the solar system (e.g. Weller & 
Meier 1981; Witte et al. 1996; Flynn et al. 1998; Frisch et al. 1999) when projected towards 
a Cen. Optical and ultraviolet (UV) data show that the bulk flow of the closest interstellar 
material (ISM) has an upstream direction towards the Loop I supernova remnant, and a 
velocity ~20 km s _1 in the local standard of rest (LSR) (Frisch 1981; Crutcher 1982; Frisch 
1995). This bulk flow is similar to ISM expanding around OB associations (e.g. Munch 1957). 
Following Slavin & Frisch (2002, hereafter SF02), the ISM within 30 pc is referred to as the 
complex of local interstellar clouds (CLIC), while the cloudlet feeding ISM into the solar 
system is denoted the Local Interstellar Cloud (LIC). Optical Ca + and UV absorption data 
show that the CLIC is inhomogeneous on subparsec scales, and that multiple absorption 
components are present towards the nearest stars (e.g. Munch & Unsold 1962; Ferlet et al. 
1986; Lallement et al. 1986; Vallerga et al. 1993; Lallement et al. 1994; Crawford & Dunkin 
1995; Welty et al. 1996; Frisch 1995, 1996; Crawford et al. 1997, 1998). The kinematics of 
the CLIC provides the opportunity to probe the history of a diffuse cloud, the relevance of 
small scale structure to ISM physics, and to gauge the past and future galactic environment 
of the Sun (Frisch 1997). 

Nearby ISM provides a unique set of constraints for determining diffuse cloud physics. 
Observations of both pickup ions inside the solar system (formed by interactions of interstel- 
lar neutrals with the solar wind) and interstellar absorption lines towards nearby stars have 
been used to constrain the first full radiative transfer model of nearby ISM (SF02). These re- 
sults indicate that the interstellar properties at the solar location are T~7,000 K, n(H°)~0.24 
cm" 3 , n(e~)~0.13 cm -3 , and fractional ionizations X(H)~31% and X(He)~48%, with both 
density and ionization levels varying towards the cloud surface (Model 17). The model in- 
cludes emission from a conductive interface, which yields an excess of helium compared to 
hydrogen ionization (SF02). If the density n(H°)~0.24 cm" 3 is typical for all cloudlets near 
the Sun, then iV(H ) towards a Cen (1.3 pc, Linsky & Wood 1996) suggests a filling factor 
/ ~0.4. Enhanced refractory element abundances in warm nearby gas provide evidence that 
local ISM has been shocked (Frisch et al. 1999). The presence of cold Ca + and Na° absorp- 
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tion components (Doppler width 6 D <0.8 km s _1 ) towards a Pav (56 pc) and 5 Cyg (52 pc) 
indicate at least an order of magnitude variation for ISM temperatures within 50 pc (see 
references in Table 3). 

The sensitivity of heliosphere properties (Zank & Frisch 1999) — and astrospheres in 
general (including extra-solar planetary systems, Frisch 1993) — to the physical conditions 
of the surrounding ISM motivates this exploration of small scale structure in the nearby 
ISM using cloud velocity as a structure proxy. Understanding small scale structure in the 
CLIC will also help decipher the signature that astrospheres leave on the interstellar H° Lya 
absorption line. This paper focuses on the bulk flow kinematics of ISM within ~30 pc of the 
Sun. Anticipating a conclusion, the CLIC appears to consist of an ensemble of cloudlets at 
velocities consistent with a random distribution about a mean bulk flow velocity. The first 
results of this study were presented by Frisch (2001). 

2. Velocities of Absorption Components 
2.1. Method 

In this analysis we assume that the motions of nearby interstellar clouds (i.e. the CLIC) 
in the local standard of rest (LSR) can be described by a linear flow vector, Vftow, which is 
characterized by the flow velocity and the galactic coordinates of the direction from which 
the gas flows (implying that the flow velocity is <0). Our approach is to calculate a best 
fitting flow vector for a set of observed interstellar cloud component radial velocities which 
sample interstellar gas within ~30 pc (Section 2.2), and then to determine the distribution 
of component velocities about this bulk flow vector. If the CLIC (d<30 pc) is part of an 
expanding shell feature from Loop I (Section 3.1), it will subtend a total angle equivalent 
to about 10% of the expanding shell and yield <5% deviations from a linear flow velocity 
(i.e. ~1 km s _1 for uniform expansion). Therefore, over the scale length of the CLIC, the 
linear assumption is sufficient. The bulk flow velocity vector is determined from a fitting 
procedure that varies Vr ow in order to minimize the sum, $ m , over m observed interstellar 
absorption line components, of the square of the difference between the projected flow and 
observed velocity towards each star, i.e., to minimize: 

m 
i=l 



where 



dVi — i>i i0bs — Vfj ow • k star . 



(2) 
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For a flow vector calculated from a set of m components (Vfl ow (m)), 

dVi(m) = v ifihs - 14ow(m) • A; star , (3) 

where dV|(m) is the scalar radial velocity (at the solar location) of component i in the rest 
frame of Vfl ow (m). Here A; star is the unit position vector towards each star and i>i )0 bs is the 
observed heliocentric (HC) radial velocity for an interstellar velocity component towards that 
star (here denoted i, where i=l,m). The fits were performed using the FindMin function 
in Mathematica, which determined the local minima of the function <3> m as a function of 
the three variables V,l,b which define Vq ow . Note that since many stars show more than 
one interstellar absorption component, the number of components (m) is greater than the 
number of stars. An unweighted fit was employed since there is no obvious weighting function. 
Line broadening is not suitable since instrument resolutions and thermal broadening vary 
(artificially blending velocity components). Column densities are unsuitable since they vary 
between species. Note that since we are interested in the distribution of the deviation of 
component velocities in the rest frame of Vfl ow , the component data set includes all measured 
components for each star. We find that this procedure yields a predicted flow velocity, Vr ow , 
which is relatively insensitive to the detailed set of the selected velocity components, provided 
that sky coverage is good and that several stars with components at anomalous velocities 
are removed from the sample (next section). 

In principle, the full three-dimensional flow vector Vfj ow is calculated correctly for cases 
where the perturbations from the flow (dVi) are random, since the Sun is immersed in the 
flow. Alternative cases whereby perturbations in the flow have a directional dependence 
(or preference), or where the flow is decelerated (Frisch 1995), are not considered here. 
However, we make an elementary attempt to identify systematic regional patterns in the 
velocity components. 

2.2. Component Set and Bulk Velocity 

The velocities of optical and ultraviolet absorption lines towards nearby stars, and to- 
wards more distant stars which sample primarily the CLIC gas, are used to evaluate the 
velocity distribution of the CLIC. Interstellar absorption line data for 67 stars which sample 
the CLIC are listed in Table 2, with component information in Table 3. This component set 
forms the "unrestricted sample" . A range of data sources was used in order to cover as much 
of the sky as possible. Observations of the optical interstellar Ca + lines represent the primary 
source of cloud velocity data, but the optical data do not adequately sample the downstream 
region where column densities are low over the first ~30 pc (e.g. iV(Ca + )<< 10 10 cm~ 2 , 
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Bruhweiler & Kondo 1982; Frisch & York 1983; Frisch 1995; Lallement et al. 1995). Ob- 
servations of the intrinsically stronger UV lines of Fe + , Mg + , D°, and H° are available for 
several stars in the downwind direction and towards the north galactic pole (where column 
densities are also low). 

The component sample results primarily from fits to optical and UV absorption lines 
observed at relatively high resolution (<3 km s _1 ); medium resolution UV data were used 
for a few stars. When both optical and UV data exist, the optical Ca + data (resolution 
0.5-3 km s _1 ) are the first choice, and the UV data second choice since about 50% of the 
optical spectra were acquired with resolution FWHM<1.5 km s^ 1 . Component velocities are 
generally obtained from a process which iteratively fits each discernible component in the 
absorption line profile with a population of atoms with a Maxwellian velocity distribution 
about the central component velocity, with a projection in the radial line-of-sight towards 
the Sun of t> bs (denoted i>i, bs in e QS. 2, 3). The number of components and the component 
descriptors (column density, t> bs, an d Doppler broadening constant, &d) are sensitive to 
instrument resolution and signal-to- noise and to the judgment of the observer. 

The component sample used here is not of uniform quality. For example, the star a 
Oph (14 pc) has been observed at resolutions of 0.3 km s _1 by Crawford & Dunkin (1995) 
and 1.2 km s -1 Welty et al. (1996). The 0.3 km s _1 resolution profiles were fitted with four 
components with velocities -32.0±0.5, -28.4±0.3, -26.2±0.2 and -23.6±0.5 km s _1 . The 1.2 
km s _1 resolution profiles were fitted with three components with velocities -33.03 km s _1 , 
-26.25 km s -1 , -22.16 km s -1 . While temporal variations in interstellar Ca + absorption pro- 
files are seen in a few cases (e.g. Hobbs et al. 1991), comparisons between the two absorption 
profiles suggest that the component at -28.4 km s -1 was unresolved in the lower resolution 
data. 

The derived flow vector depends on the component set used in the fitting procedure, 
and the unrestricted sample is small enough that lines incorrectly identified as interstellar 
may alter the results. Therefore, the basic star sample was restricted to omit stars with 
either identified circumstellar disks or strong emission line spectra (where misidentification 
of circumstellar features is possible). Thus the stars e Eri (3.2 pc), (3 Pic (19.3 pc), (3 Leo 
(11.0 pc), a Lyr (7.6 pc), a PsA (7.7 pc), which have circumstellar disks identified at 60 /im 
(e.g. Habing et al. 2001), were omitted from the sample. Although interstellar components 
can in principle be distinguished from circumstellar disk components since stellar velocities 
are known, we considered the sample to be more reliable when stars with circumstellar disks 
are omitted. The emission line star a Eri (44 pc) was also omitted, as it shows variable He° 
and Mg + chromospheric features indicating that weak circumstellar absorption features may 
be present. The active RS CVn variable HR 1099 (Wood et al. 1996b) was also omitted from 
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the sample because the inclusion of the components towards this star yielded results that 
are inconsistent with fits towards the other stars (see below). With the above omissions, 
the remaining restricted sample consists of 60 stars with a total of 96 components. Star 
distances range from 1.3 pc (a Cen) to 132 pc (e CMa). Components towards the two stars 
beyond 70 pc, e CMa and 31 Com, were included because these stars primarily sample ISM 
within 30 pc (Gry & Jenkins 2001; Dring et al. 1997) and fill a gap in the spatial coverage. 

This basic component sample was fitted for the flow vector Vfl OW) yielding a HC vec- 
tor Vfl ow (96) with radial velocity -28.1±4.6 km s _1 , flowing from the upstream direction 
/ n =12.4°, 6 n =11.6°. The uncertainty is the 1-a value of the velocity component distribution 
about the bulk flow velocity. Figure la shows the deviations dVi(96) plotted against the 
stellar distance (omitting the two stars more distant than 70 pc). Note that 76 of these 96 
components originate in the upstream direction, while the remaining components originate 
downstream, reflecting the higher column densities of nearby ISM found upstream (e.g., 
Frisch 1995). 

The positions of the stars forming the 96 component set are plotted in Figure 2a. Most 
of the sky is well sampled, except for the interval Z n =150°-180 ° in the direction of the third 
galactic quadrant void (also known as the Puppis window or the (3 CMa tunnel), where 
mean interstellar column densities are low out to ~100 pc (njjo <0.003 cm -3 ). Figure 3a 
shows dVi(96) as a function of the galactic longitude of the star; no strong systematic spatial 
dependence is seen. For comparison, the LSR velocities of the 96 components (dVj(LSR)) 
are plotted against galactic longitude in Figure 3b, illustrating the expected behavior for an 
incorrect vector describing the flow. Throughout this paper we use the solar apex motion 
based on Hipparcos data, corresponding to a solar motion towards the galactic coordinates 
/ n =27.7°, 6 n =32.4° at the velocity V=13A km s" 1 (Dehnen & Binney 1998). 

The robustness of Vfj ow (96) was tested in several ways. The first test was to restrict the 
star sample to stars within 29 pc of the Sun (except for e CMa and 31 Com, which sample 
the nearest ISM and are required for completeness in low column density directions). This 
restriction gives 46 components towards 31 stars. The best fit HC flow vector for this 46 
component sample, Vq ow (AQ) : corresponds to an upstream direction of / n =12.5°, 6 n =12.5°, 
with an inflow velocity of-28.1±4.3 km s _1 . The fact that Va ow (46) is virtually identical to 
Vfiow(96) is not a coincidence, since the initial selection of stars in Table 3 was restricted to 
stars which did not appear to contain neutral ISM more distant than ~30 pc. The closeness 
of the l / How(96) and Vfl ow (46) vectors support the validity of eliminating HR 1099 from the 
component sample, since had HR 1099 been retained in the restricted sample the resulting 
fitted vector (Vfl ow (99)) would have differed from \? fl ow(46) by 0.6 km s _1 and 2°. 

The second test was a search for spatially correlated variations in dV;(96). Star posi- 
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tions are replotted in Figure 2b, with symbols coded for dVJ(96). No systematic positional 
dependence appears for components with | dV^ (96) | >5 km s -1 . However, when the compo- 
nents with |dV;(96)| <5 km s _1 are plotted with symbols which code the sign of dVi(96), 
stars in the interval / n =150°^250° are found to show components with systematically small 
but negative dVi(96) values (Figure 2c, note several components are overp lotted) . The best- 
fitting flow vector, Vfl ow (20), for the 20 components toward 14 stars located in the interval 
/ n =150 o -^250° corresponds to a vector velocity -25.8±4.3 km s _1 from the upstream direc- 
tion / n =6.2°, 6 n =10.4°. If only a single component per star is selected for the fit, biasing 
towards components with small d^(He) values, the resulting fitted vector V^ ow (14) cor- 
responds to a velocity -25.6±1.3 km s^ 1 from the upstream direction / n =6.4°, 6 n =13.1°, 

— * 

which is close to Vfl ow (He) (Table 1, Section 3.3.1). The dispersion of the same 14 compo- 
nents around \4 ow (He) is ±1.7 km s -1 . Gas near the LIC velocity (Table 1) dominates in the 
downwind direction (1~180°), as found previously. The LIC and other regional variations in 
Vaow(96) are discussed in Section 3.3. 



2.3. Component Velocity Dispersion 

The randomness of interstellar cloud velocities has long been of considerable interest, 
with early arguments presented for random motions of ~20 km s -1 for interstellar particles 
(Spitzer 1942), compared to recent data showing root-mean-square velocity dispersions of 
~0.5 km s _1 in individual cold clouds. Early studies of diffuse cloud kinematics found an 
exponential distribution for cloud velocities, suggestive of a turbulent ISM (Blaauw 1952; 
Munch 1957). 

The functional form of the velocity distribution of the CLIC gas was tested, but the 
results prove inconclusive. A plot of the number of components (N, ordinate) binned for a 
given value of dVJ(96) (abscissa) is shown in Figure 4. Here, N is determined with dV^(96) 
binned with 1 km s -1 bin sizes. The histogram is "noisy" because the sample is small. 

For a purely random distribution of individual components about the central flow ve- 
locity Vfl ow (96), the form of the dVi distribution should be Gaussian: 

= T^exp - (V - V f/tf (4) 

The Blaauw (1952) and Munch (1957) results for Ca + and Na° lines, based on low resolution 
data (FWHM>5 km s _1 ), are consistent with an exponential distribution for cloud velocities: 

V(y) = £-exp-\V-V \/ri, (5) 
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such as is expected for a turbulent flow. Here, b = c g and c c are normalizing con- 

stants, and a and r)\/2 are the root-mean-square deviations of the Gaussian and exponential 
distributions respectively. 

The 96 component sample was tested for the distribution of dV|(96), assuming alter- 
natively Gaussian and exponential distributions. The best-fit Gaussian distribution yields 
6 n =6.2 km s _1 , and the best-fit exponential distribution yields 77=5.1 km s _1 . The present 
data do not distinguish between these two distributions (see overplotted functions in Fig- 
ure 4). Blaauw concluded that an exponential function with mean speed rj ~5±1 km s _1 
provided the best fit to Ca + K-line data towards 300 stars observed by Adams (1949) with 
instrumental resolution ~9 km s _1 . Munch inferred that an exponential form fit observations 
of the Ca + and Na° doublets towards 112 stars, with rj ~5-6 km s _1 for Ca + in the Orion 
arm and a somewhat smaller value for Na°. Munch concluded that turbulence dominates 
the velocity distribution of these components. These early low resolution data undersample 
blended components at low velocity and therefore probably enhance component statistics in 
the distribution wings with respect to the central region. It is therefore puzzling that we find 
a similar value from higher resolution data. This issue is discussed further in a subsequent 
paper (Frisch & Welty 2002). 



3. Discussion 

3.1. Upstream Direction and Loop I 

The best-fit heliocentric velocity vector, Vfl ow (96), corresponds to a bulk velocity -28.1 
km s _1 flowing from the upstream direction / n =12.4°, 6 n =11.6° (Section 2.2). However in 
order to compare Vfl ow (96) with morphological features in the astronomical sky the solar 
apex motion through the LSR must be subtracted. Subtracting the solar apex motion from 
Vfk>w(96) gives a "true" upstream direction (in LSR) in galactic coordinates of Z n =2.3°, 
b n =-5.2°, with flow velocity -17.0 km s -1 . 1 The Vfl ow (96) upwind direction is shown plotted 
against the filamentary H° structure associated with Loop I (21 cm data, Figure 5, Hartmann 
& Burton 1997). The position of the radio continuum Loop I shell, defined from the discovery 
408 MHz radio continuum data, is marked (from Berkhuijsen 1971; Haslam et al. 1982). 

Low column densities in the CLIC components (iV(H)< 10 18 ' 5 cm -2 , typically) prevent 



1 For comparison, the standard solar apex motion (based on the brightest stars with a range of ages) is 
19.7 km s _1 towards Z n =57°, 6 n =+22°, which yields an upstream LSR flow motion of -19.4 km s -1 from 
the direction Z n =331.4°, 6 n =-4.9°. 
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linking CLIC velocities to individual features in the 21-cm H° sky maps. The relation between 
the flow of interstellar gas past the Sun and the Loop I superbubble has been discussed 
from several viewpoints (e.g. Frisch 1979, 1981; Crutcher 1982; Frisch 1995; Crawford 1991; 
de Geus 1992). The Loop I superbubble remnant (the brightest segment of which is the 
North Polar Spur) is centered near Z n ~320°, 6 n ~5°, with distance ~ 130 pc and radius 
~80° (based on H° 21 cm data, Heiles 1998). The radio-continuum Loop I is centered at 
Z n =329°, 6 n =+17.5°, radius 58°, and is strongly limb-brightened in the tangential direction 
near Z n ~35° (Berkhuijsen 1973). The center of Loop I defined from H° data is offset from 
the radio-continuum center by ~15°. The H° shell shows a polar cap at t>i sr ~-30 km s" 1 and 
F,6 n a;300 o ,-10 o (Heiles 1989), where the cold component seen at -19.6 km s -1 towards a 
Pav may originate (Table 3). This expansion velocity implies a projected velocity of «-14 
km s _1 for the Vflow(96) LSR upwind direction, versus the best-fit LSR velocity of -17.0 
km s -1 . The filamentary structure and likely asymmetric expansion of the H° shell make it 
difficult to identify the approaching portion of the H° shell, so the ~20% difference between 
the projected polar cap velocity and the LSR bulk velocity of the CLIC gas possibly is 
within uncertainties. However, the H° shell regions traced by the 21-cm polar cap emission 
are denser and colder than the CLIC gas appears to be. Alternatively, models of shell 
expansion from the formation of the Sco-Cen Association Upper Scorpius subgroup place an 
H° shell with LSR velocity ~-22 km s -1 and age ~4 Myrs at the solar location (Crawford 
1991; de Geus 1992; Frisch 1995). 2 



3.2. Comparison with Previous Models 

The bulk flow of nearby interstellar gas has been determined from optical absorption 
lines observed towards nearby stars (primarily Ca + , e.g. Crutcher 1982; Frisch 1986; Lalle- 
ment et al. 1986; Lallement & Bertin 1992; Vallerga et al. 1993; Frisch 1995, 1997). The 
velocity vectors determined from these earlier studies are summarized in Table 5 and shown 
in Figure 7. Since column densities in the downstream direction are low (log[iV(H°)]<18 
cm -2 ), optical lines do not easily sample this interval. This causes bulk flow vectors based 
primarily on Ca + data to be strongly weighted towards upstream gas and quite sensitive to 
the details of the star sample. The inclusion of UV data provides an adequate sample of 
downstream gas, but at lower resolution (3 km s -1 ). Lallement & Bertin (1992, LB92) used 
UV data to identify the ISM in the downstream direction, denoting it the "anti-Galactic" 
cloud (AG, Table 5), which is the same as the LIC. ISM identified in the upstream direction, 



2 Note the age given for the most recent shell-forming event should be 4 Myrs, not 400,000 yrs (Frisch 
1995, page 532). 
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generally towards the galactic-center hemisphere, was denoted the "Galactic" cloud (G, Ta- 
ble 5). The velocity difference between the G and AG clouds is smaller than found by Adams 
and Frisch (1977), based on backscattered La H° data, because of the deceleration and com- 
pression of H° in the heliosheath regions (which was unknown in 1977). LB92 concluded that 
either the Sun is located in a patch of gas with velocity intermediate between the G and AG 
clouds, or that H° is decelerated in the heliosheath (or both). The Ulysses He and 584 A 
backscattering data confirm the velocity of ISM inside the solar system, so that the general 
velocity difference between the upstream and downstream ISM is real. However, the fact 
that the velocities of CLIC cloudlets are consistent with a random distribution about the 
mean bulk flow velocity suggests that the upstream/downstream velocity gradient locally 
displayed by the G versus LIC vectors may be an artifact of the location of the Sun with 
respect to the "edge" of the CLIC complex. Alternatively, the presence of a blue-shifted 
cloud towards CMa (Lallement et al. 1994; Gry & Jenkins 2001) is also consistent with a 
velocity gradient. 

For comparison with the LB92 G-AG two-flow model, a separate bulk flow velocity was 
calculated for downstream components. The 20 LIC components originating in the inter- 
val / n =150 o ^250° were removed from the 96 component sample. The fit to the resulting 
76 components yields \4ow(76)=-29.3±4.0 km s _1 , from HC upstream direction / n =13.1°, 
6 n =11.2°. The velocities of Vr w(7Q) and the G vector (Vfl ow (G)) coincide, although the 

— * 

upstream directions differ by ~12°. The observed 3.5 km s difference between Vfl ow (76) 
and V Row (20) and the near coincidence of upwind directions (within ~8°) suggests a decel- 
eration of the bulk flow. The two-flow model with the two vectors V aow (20) and Vfl ow (76) 

— * — * 

(or alternatively with Vfl ow (He) and Vfl ow (G)) is not a unique description of CLIC gas. For 
example, the velocity dispersion of Vfl ow (76) is about equal to the difference between the 
^4ow(76) and \4 OW (20) velocities (~4 km s _1 ). The second shortcoming is that 6 stars in 
the / n =150 o ^250° interval show two absorption components, and the second components 
are not accommodated by the two-flow model. The dispersion of the entire 20 components 
about Vfl ow (96) is 5.4 km s -1 . The dispersion around U flow (He) of the 14 components best 
matching the LIC velocity (1.5 km s" 1 ) is only slightly better than the dispersion of these 
components about Vfl ow (96) (1.6 km s^ 1 ). 



3.3. Regional Properties 

The small filling factor of nearby ISM, combined with the component distribution, led to 
attempts to understand the positions of nearby 'cloudlets' (e.g. Lallement et al. 1986; Frisch 
1996), the clouds towards a Oph and a Aql (Munch & Unsold 1962; Frisch 1981; Ferlet et al. 
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1986), and the "shape" of the cloud around the solar system (Frisch 1996; Redfield & Linsky 
2000). These earlier approaches used velocities grouped either in the observed heliocentric 
velocity, or in the LIC velocity. Here we use an alternative approach, and evaluate cloud 
membership in the rest frame of Vfl ow (96). Using the dV|(96) values, we identify regional 
clumps of ISM based on similar dVj(96) values for stars sampling a relatively compact region 
of the sky, except for the LIC which is identified by dVJ(He) (or dV|(20)). The component 
groups that appear to represent subsets of the flow are listed in Table 4 and are discussed in 
order below. The positions of the clouds are plotted in Figure 6 as symbols overlying star 
positions. 



3.3.1. 1. Cloud Surrounding the Solar System 

The velocity of the interstellar cloud surrounding the solar system, the LIC, has been 
determined from observations of the resonance fluorescence of the 584 A transition from 
interstellar He° inside the solar system (e.g. Weller & Meier 1981; Flynn et al. 1998) and from 
Ulysses in situ measurements of interstellar He° (Witte et al. 1996). Since the trajectories 
of interstellar H atoms inside the solar system are subject to radiation pressure, these He° 
data yield the best LIC velocity. The Ulysses GAS-detector data give Vfl ow (He)=-25.3±0.5 
km s _1 , Z n =3.9±1.0°, 6 n =15.8±1.3° (and T=6,550±l,050 K) for the HC upstream direction 
(Witte et al. 1996, and private communication), which is within the uncertainties of the 
He° 584 A backscattered radiation vector (Flynn et al. 1998). When this He motion is 
converted to a vector in the local standard of rest (LSR) using the Hipparcos solar apex 
motion (Dehnen & Binney 1998), the upstream direction corresponds to a velocity of -14.7 
km s _1 , arriving from the upwind direction / n =345.9°, fe n =-l.l°. In the LSR, this differs 
by ~2 km s _1 and ~20° in direction from the bulk velocity vector. Pure deceleration of the 
flow would not alter the upstream direction, so shear motions may be present in the bulk 
flow. 

Copernicus, IUE and HubbleSpaceTelescope (HST) observations of interstellar D°, 
Fe + , Mg + , towards a Cen (1.3 pc) indicate a HC velocity for the ISM in the range -17.7±0.1 
to -18.2±0.1 km s _1 , versus the value -16.1 km s _1 predicted by the projection of the 
Ulysses vector (Landsman et al. 1984; Lallement et al. 1995; Linsky & Wood 1996). This 
difference has been interpreted as indicating that the cloud surrounding the solar system 
terminates within ~10,000 au of the Sun in the direction of a Cen. The a Cen cloud shows 
dVi(He)~-2.8 km s _1 (and dVJ(96)=-3.8 km s _1 ). This evidence for a cloud boundary close 
to the Sun in the direction of a Cen is consistent with column densities measured towards 
the nearest stars, which indicate that <40% of space is filled with ISM for volume density 
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n(H°)~0.24 cm" 3 as found for the LIC (SL02). 

Thirty CLIC components have velocities consistent with the velocity of interstellar He 
inside the solar system (V flow (He)) using the criterion |dVj(He)| <1.5 km s _1 . However, 
components towards stars beyond 15 pc (Figure lb) have predominantly negative dVi(He) 
values since these stars are generally located in the upstream direction. The stars showing 
these 30 Vfl ow (He) components are distributed relatively uniformly across the sky (Figure 
2d). Refitting the 30 components essentially reproduces Vfl ow (He), and does not alter the 
distance dependence, suggesting the more distant components with |dVi(He)| <1.5 km s -1 
either are blends of local and distant gas or are formed in independent cloudlets with no local 
contribution. The five stars within 15 pc with Va ow (He) components are a CMa (2.7 pc), 61 
CygA (3.5 pc), e Ind (3.6 pc), 40 EriA (5 pc), and a Aur (13 pc). The UV absorption lines 
towards a CMi (3.5 pc) have been fitted with both single- component and two-component 
models (Linsky et al. 1995). With the single-component fit (Table 3), dVi(He)=1.8 km s _1 ; 
however dVJ(20)=0.5 km s" 1 so that Vfl ow (20) provides a better fit for local gas towards this 
star. (A second component redshifted by ~2.6 km s _1 (with ~50% of the column density of 
the main component) is not coincident with either flow vector.) 

The velocity vector of the cloud within the solar system provides poor agreement with 
the restricted sample of 96 components, as seen in Figure 3c where dVj(He) is plotted for all 
96 components. The fact that |Va ow (He)| < |Vfl ow (96)| yields the effect that nearby ISM in 
the downstream direction is blueshifted in the rest frame of Vfl ow (96). 

Several studies have noted the small distance to the upstream cloud edge for the ISM 
surrounding the solar system (< 10 4 au, e.g., LB92). The relative Sun-cloud velocity (~25 
pc/10 6 yrs) suggests that the galactic environment of the Sun will vary within the next 
~2,000 years (Frisch 1997; Wood et al. 2000b). 

The LIC is reliably identified only in the / n =150°— >250° interval where slightly more 
distant gas is absent, and in this interval Va w(He) is an excellent descriptor of component 
velocities. The inability to identify the LIC reliably in sidestream and upstream directions 
indicates the two-flow model must be evaluated with caution in most sightlines, and in 
particular where velocity information is unavailable. In addition, since the CLIC represents 
a cluster of comoving clouds with a velocity dispersion of ~5 km s -1 , simple closed surface 
LIC models which can be topologically deformed to a sphere may not realistically represent 
the LIC cloudlet. 
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3.3.2. Other Regional Component Groups 

2. Blue Cloud: The "Blue Cloud" (BC) has been identified towards a CMa and e CMa, 
placing the cloud within 3 pc of the Sun (Lallement et al. 1994; Hebrard et al. 1999; Gry 
& Jenkins 2001, GJ). The BC has dV;(96)~-10 km s _1 , and is blueshifted from the LIC by 
~6.5 km s _1 . The filling factor of the LIC towards a CMa is / ~ 0.3, suggesting that the 
BC is a spatially separate feature. Towards a CMa the BC appears to be cooler and denser 
than the LIC (Hebrard et al. 1999), while towards e CMa it appears warmer and denser than 
the LIC (GJ). 

3. Aquila-Ophiuchus Cloud: The Aquila-Ophiuchus cloud (AOC) is disclosed by a 
set of Ca + components with dV^(96)=-6.1±0.9 km s _1 towards stars located in the interval 
/ n =28° to 48°, 6 n =-5° to +23°. The stars a Aql, a Oph, ( Aql, 7 Oph, and A Aql show 
the Aql-Oph cloud component, and the 5.1 pc distance of a Aql places the cloud close to 
the Sun. For each of these five background stars, the component formed in the AOC is the 
weakest and most blue-shifted component seen towards the star. The similarity between the 
AOC velocity and the velocity of weakest and most blueshifted component towards 5 Cyg 
(dVi(96)=-6.9 km s" 1 ) may be a coincidence since the the 5 Cyg component originates in a 
cold cloud (6=0.47 km s -1 , 0.42 km s" 1 for Ca+, Na° respectively, Welty et al. 1996), and 
cold clouds are not expected within 5 pc. The AOC is located near the solar apex direction, 
/ n =27.7°, 6 n =32.4°. If the unobserved tangential velocity component is ~0 km s -1 , the Sun 
may encounter this cloud within the next ~160,000 years. Components from this cloud are 
seen clearly in Figure 3a. 

4. Pegasus-Aquarius Cloud: The Pegasus-Aquarius cloud, with dV[(96)=4.2±0.5 
km s _1 , is seen towards stars grouped in the interval / n =75°±13 , 6 n =-44°±5°. Background 
stars are 9 Peg, a Peg, 7 Aqr, 77 Aqr, and this cloud must be within 30 pc of the Sun. 

5. North Pole Cloud: A group of components with velocity 1.7±0.6 km s~ x is found 
towards stars at high galactic latitudes, 6 n >53°, including the Ursa Majoris region. The 
nearest star with a component in this group is a CrB, at 23 pc. 

6. South Pole Cloud: A cloud with \? fl ow(96)=2.4±1.0 km s _1 occupies an irregularly 
shaped region covering the South Galactic Pole (6 n <48°). The South Pole Cloud (SPC) cloud 
is seen towards e Ind, a Hyi, r 3 Eri, (3 Cet. However, two of these components (towards e 
Ind and (5 Cet) also have |dV;(He)| <0.1 km s -1 , indicating they also could be formed in the 
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LIC (which is seen towards two other low latitude stars, 40 Eri and EP Eri). If so, the SPC 
is not a separate cloudlet. The nearest star in this sample is e Ind at 3.6 pc, indicating the 
SPC, if real, is quite nearby. Additional data are required to test for the SPC. 

7. Filamentary-like Feature: Perhaps the most puzzling of the component groups 
suggests a filamentary structure, with dVj(96)= +6.9±1.5 km s _1 . Since this feature has a 
large angular extent (~90°, Figure 6), it represents a single cloudlet only if the tangential 
(unobserved) velocity is ~0 km s -1 for all of the stars, similar to a fragment of an expanding 
ring. This feature is at Z n =284°±12°, 6 n =-4°±50°, and it is seen towards 5 CrV, 7 CrV, 
HR 4023, 5 Vel, a Hyi. 

4. Conclusions 

The principal conclusions of this paper are that: 

1. The bulk flow velocity for interstellar matter within ~30 pc of the Sun is determined 
from optical and UV absorption line data. A self-consistent flow vector is determined if 
stars with infrared-emitting circumstellar disks (e Eri, (3 Pic, (3 Leo, a Lyr, and a PsA), 
variable emission lines (a Eri), and HR 1099 are omitted from the sample. Fits to the 
remaining 96 component sample (60 stars) yield a bulk flow velocity (Vfl ow (96)) through the 
LSR (using Hipparcos data on the solar apex motion) of magnitude -17.0 km s" 1 , with an 
upstream direction of Z n =2.3°, 5 n =-5.2° (Table 1). In the heliocentric velocity frame (i.e. 
with respect to the Sun), |Vfl ow (96)|=28.1 km s _1 , approaching the Sun from the upstream 
direction Z n =12.4°, ft n =+11.6° (the l-a velocity uncertainty is ±4.6 km s _1 , based on the 
component distribution about the bulk flow velocity). 

2. The CLIC gas (Cluster of Local Interstellar Clouds) is defined by this ensemble of 
velocity components, which show a Gaussian or exponential distribution in the rest frame 
of the central flow velocity, and the dispersion of this distribution is typical for dispersions 
determined from lower resolution observations of diffuse clouds. The dispersion in the veloc- 
ities of the individual CLIC components partaking in the bulk flow explains the difference 
between the velocity of the interstellar cloud inside the solar system (25.3 km s -1 , HC) versus 
the overall bulk flow velocity (-28.1 km s _1 , HC). 

3. This LSR upstream direction suggests the CLIC gas may be part of a superbubble 
shell expanding from the Loop I supernova remnant, or from an earlier epoch of superbubble 
formation in the Scorpius-Centaurus Association. The cold absorption component formed 
in the polar cap of Loop I appears to be present in front of a Pav. 



- 15 - 



4. The velocity of the interstellar cloud observed inside the solar system, and in the 
downstream direction towards nearby stars (Vflow(He)) differs both from Vfl ow (96), and from 
the velocity of ISM towards the nearest star a Cen. However, 85% of the observed compo- 
nents consistent with Vfj ow (He) are probably formed in unrelated parts of the CLIC. This 
complexity suggests that simple smooth closed-surface models for local interstellar gas are 
unlikely to be accurate in detail. The components within 1.5 km s _1 of the velocity of 
interstellar He within the solar system (\? fl ow(He)) dominate the galactic longitude interval 
/ n =150°— ^250°, suggesting a deceleration of the flow in the third galactic quadrant. The Sun 
is likely to emerge from the gas-patch now surrounding it within the next several thousand 
years. 

5. Several spatially distinct groups of components sharing a common dV|(96) suggest the 
bulk flow is composed of cloudlets (Table 4), including the previously known LIC and Blue 
Cloud in the downstream direction. Distinct cloudlets in the upstream direction include the 
Aquila-Ophiuchus cloud with components at at dVj(96)~-6 km s -1 . This cloud may extend 
to within 5 pc of the Sun since it is seen towards a Aql. Components in this cloud form 
the weakest and highest velocity Ca + component for each star in which it is observed. This 
cloud is in the solar apex direction, and if the non-radial component of the cloud velocity is 
small the Sun will encounter this cloud in <200,000 years. A distinct cloud is seen towards 
Pegasus/Aquarius, at dV;(96)=3.8±0.3 km s" 1 . Another distinct cloud is found within 22 
pc towards the North Galactic Pole, with dV^(96) = 1.7±0.6 km s -1 . 

6. These results show the kinematics of CLIC components are consistent with early 
studies showing macroscopic turbulence in the ISM within ~500 pc (Blaauw 1952; Munch 
1957). Low filling factors for nearby ISM(<40%), and mean cloudlet speed rj ~5 km s _1 in 
the flow reference frame (t4 ow (96)) suggest that the CLIC is a cluster of comoving cloudlets. 
If this velocity dispersion is due to turbulence (as the larger samples considered by Blaauw 
and Munch suggest), then it is distinct from the microscopic turbulence which broadens 
absorption lines towards nearby stars (b 2 =2kT/m + 2V^ ulh , where \4urb ~1 km s _1 ). 

7. A two-flow model of the velocity components is not a unique description of CLIC 
gas. For example, the velocity dispersion of Vfl ow (76) is about equal to the difference between 
the Vfk>w(76) and Vfl ow (20) velocities (~4 km s _1 ). In addition, multiple components in the 
downstream gas are not accommodated by a two-flow model. The result is that morphological 
models of the LIC must include velocity information for accurate results. 

Studies of the CLIC, where individual cloudlets can be studied in high resolution UV 
data, and in some sightlines in optical data, offer the best opportunity for reliably deter- 
mining the relation between the kinematical and spatial characteristics of diffuse interstellar 
clouds. 
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Table 1. Streaming Motions in HC and LSR Rest Frames 



HC Vector a LSR Vector a ' b Sample c 



Vel. 1 b Vel 1 b 

km s _1 Deg. Deg. km s _1 Deg. Deg. 



Bulk Fits: 
















Fflow(96) 


-28.1±4.6 


12.4 


11.6 


-17.0 


2.3 


-5.2 


60 stars 


F fl ow(46) 


-28.1±4.3 


12.5 


12.5 








31 stars, d<29 pc 


Upstream: 
















F flow (76) 


-29.3±4.0 


13.1 


11.2 


-18.1 


3.9 


-4.7 


46 stars 


Fflow(G) 


-29.4 


4.5 


20.5 








G-cloud (Lallement & Bertin 1992) 


Downstream, 


Z n =150^250 














F flow (He) 


-25.3±0.5 


3.9±1.0 


15.8±1.3 


-14.7 


345.9 


-1.1 


Ulysses He° value 


F flow (20) 














(Witte et al. 1996) 


-25.8±4.3 


6.2 


10.4 








14 stars 


F flow (584A) 


-26.4±1.5 


4.5±0.6 


14.2±0.6 








Backscattered He° 584 A 
















(Flynn et al. 1998) 



a Galactic coordinates correspond to upstream directions. The velocity uncertainty is the 1-a value of the compi 
nent distribution about the bulk flow velocity. 

b Based on Hipparcos solar apex motion (Dehnen & Binney 1998). 
c Sample stars are listed in Tables 2 and 3. 



Table 2. Data Sources 



Reference Nominal Star List 

Resl. a 
km s _1 



Optical Data: 

Frisch & Welty (2002, FW02) 1.3-1.5 

Crawford et al. (1998, CLW98) 0.35 

Crawford et al. (1997, CCW97) 0.35 

Welty et al. (1996, W96) 1.2 
Crawford & Dunkin (1995, CD95) 0.3 

Vallerga et al. (1993, V93) 1.9 

Bertin et al. (1993, B93) 3.0 

Lallement & Bertin (1992, LB92) 3.0 

Ultraviolet Data: 

Gry & Jenkins (2001, GJ01) 2.7 

Wood et al. (2000b, WLZ00) 2.7 

Wood et al. (2000a, WOO) ~15 

Sahu et al. (1999, Sa99) 2.7 

Hebrard et al. (1999, H99) 2.7 

Wood & Linsky (1998, WL98) 3.5 

Dring et al. (1997, D97) 3.5 

Wood et al. (1996a, WAL96) 3.5 

Piskunov et al. (1997, P97) 3.5 

Ferlet et al. (1995, Fer95) 3.5 

Linsky et al. (1995, Lin95) 3.5 

Lallement et al. (1995, Lal95) 3.5 



a CMa, a Leo, rj UMa, 7 UMa, 5 UMa, (3 Lib 

a CrB, A Oph, 5 Her, v Ser, (3 Ser, o Ser, A Aql 

C Aql, a Aql, (3 Aur, a Peg 

a And, 5 Vel, a Pav, 6 Peg 

a Hyi, 1 Cen, e Gru, 7 Aqr, 7 Oph 

7 CrV, S CrV, 6 Cyg 

a Oph, a Gru, a Eri 

al CVn, 7] Aqr, k And 

r3 Eri, HR 4023 

a Cep, 5 Cas, a Lac 

e CMa 
36 Oph 

EP Eri, DX Leo, LQ Hya, V368 Cep 

G191-B2B 

a CMa 

61 Cyg, 40 Eri A 

e Eri, (3 Gem, a Tri, a Gem, 31 Com 
e Ind, A And 
13 Cet 
a PsA 

a CMi, a Aur 

a Lyr, (3 Leo, (3 Pic 



-22 - 



Table 3. Velocity Components of Unrestricted Sample 



HD 


Name 


1 


b 


d 


Spec. 


N(X) a 


Velocity 


dVi(96) 


Ref. b 






deg. 


deg. 


pc 






km s _1 


km s _1 




128621 


a CenB 


315.7 


-0.7 


1.4 


Glpl 


3.89el7 


-18.9 


-3.8 


P97 


48915 


a CMa 


227.2 


-8.9 


2.6 


A1V 


2.5el7 


13.7 


-9.5 


H99 














0.16el0 


19.6 


-3.6 


FW02 


22049 ' 


' T-l ■ 

e En 


196.0 


-48.0 


3.2 


K2V 


7.50el7 


21.3 


-1.3 


D97 


61421 


a CMi 


213.7 


13.0 


3.5 


F5IV/V 


7.59el7 


20.8 


-3.0 


Lin95 


201092 


61 CygA 


82.3 


-5.8 


3.5 


K5V 


7.08el7 


-9.0 


0.0 


WL98 














7.41el7 


-3.0 


6.0 


WL98 


209100 


e Ind 


336.2 


-48.0 


3.6 


K4.5V 


1.00el8 


-9.2 


1.5 


WAL96 


26965 


40 EriA 


200.8 


-38.0 


5.0 


K1V 


8.71el7 


21.6 


-3.4 


WL98 


155886 


36 OphAB 


358.3 


6.9 


6.0 


K1V/K0V 


7.08el7 


-28.4 


-1.2 


WLZ00 


187642 


A 1 

a Aql 


47.7 


-8.9 


5.1 


A7V 


0.3el0 


-26.9 


-5.5 


FW02 














0.5el0 


-22.8 


-1.4 


FW02 














0.9el0 


-18.1 


3.3 


FW02 


172167 


d t 

' a Lyr 


67.5 


19.2 


7.6 


A0V 


6.5el2 


-18.3 


-1.5 


Lal95 














6.9el2 


-16.0 


0.8 


Lal95 














1.03el3 


-12.7 


4.1 


Lal95 


216956 c 


• a a PsA 


20.5 


-64.9 


7.7 


A3V 


3.8el2 


-3.2 


3.3 


Fer95 














1.8el2 


-9.7 


-3.2 


Fer95 


17925 


EP Eri 


192.1 


-58.3 


10.4 


K2V 


8.9el7 


19. 


-0.3 


WOO 


62509 


/3 Gem 


192.2 


23.4 


10.6 


K0II 


1.15el8 


21.9 


-1.2 


D97 














6.82el7 


33.0 


9.9 


D97 


102647 c 


' d /3 Leo 


250.6 


70.8 


11.0 


A3V 




-0.8 


-0.3 


Lal95 
















2.7 


3.2 


Lal95 
















13.3 


13.8 


Lal95 


34029 


a Aur 


162.6 


4.6 


12.9 


G5IIIepl 


1.74el8 


22.0 


-1.4 


Lin95 


159561 


a Oph 


35.9 


22.6 


14.3 


A5III 


l.OOell 


-23.6 


1.9 


CD95 














1.58ell 


-26.2 


-0.7 


CD95 














1.59el0 


-32.0 


-6.5 


CD95 














6.31el0 


-28.4 


-2.9 


CD95 


203280 


a Cep 


101.0 


9.2 


15.0 


A7IV 


0.7el0 


0.2 


1.8 


LB92 


432 


/3 Cas 


117.5 


-3.3 


16.7 


F2IV 


1.51el8 


10 


2.5 


P97 
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Table 3 — Continued 



HD 


Name 


1 


b 


d 


Spec. 




Velocity 


dVi(96) 


Ref. b 






dee. 


dee. 


DC 






km s ^ 


km s ^ 




11443 


a Tri 


138.6 


-31.4 


17.5 


F6IV 


1.15el8 


17.6 


0.8 


D97 














9.79el7 


12.7 


-4.1 


D97 


82443 


DX Leo 


201.2 


46.1 


17.7 


K2V 


5.el7 


11. 


-3.8 


WOO 


115892 


i Cen 


309.4 


25.8 


18.0 


A2V 


2.5el0 


-18.2 


-4.5 


CCW97 


82558 


LQ Hya 


244.6 


28.4 


18.3 


K2V 


5.62el8 


6. 


-6.2 


WOO 


39060 c ' d 


(3 Pic 


258.4 


-30.6 


19.3 


A5V 




-3.2 


-15.7 


Fer95 
















-9.7 


-22.2 


Fer95 


220140 


V368Cep 


118.5 


16.9 


19.7 


K2V 


8.9el7 


5. 


-0.6 


WOO 


12311 


a Hyi 


289.5 


-53.8 


21.9 


F0V 


2.88el0 


9.8 


7.3 


CCW97 














5.01el0 


4.9 


2.4 


CCW97 


139006 


a CrB 


41.9 


53.8 


22.9 


A0V 


1.03el0 


-17.4 


1.3 


FW02 


87901 


a Leo 


226.4 


48.9 


23.8 


B7V 


0.5el0 


10.5 


-0.3 


FW02 


156164 


5 Her 


46.8 


31.4 


24.1 


A3IV 


3.1el0 


-19.7 


2.7 


FW02 


74956 


5 Vel 


272.1 


-7.4 


24.4 


A1V 


0.86el0 


15.6 


10.0 


CLW98 














2.38el0 


1.3 


-4.3 


CLW98 














4.34el0 


11.6 


6.0 


CLW98 


106591 


5 UMa 


132.6 


59.4 


25.0 


A3V 


0.9el0 


3.8 


1.6 


FW02 


40183 


13 Aur 


167.5 


10.4 


25.2 


A2IVpl 


l.OelO 


22.3 


-1.3 


FW02 


177724 


C Aql 


46.9 


3.3 


25.5 


AOVn 


0.6el0 


-24.5 


-1.5 


FW02 














0.9el0 


-20.6 


2.4 


FW02 














1.4el0 


-29.7 


-6.7 


FW02 


103287 


7 UMa 


140.7 


61.4 


25.6 


A0V 


0.8el0 


4.4 


1.2 


FW02 


222107 


A And 


109.9 


-14.5 


25.8 


G8III 


2.81el8 


6.5 


1.6 


WAL96 


18978 


r3 Eri 


213.5 


-60.3 


26.4 


A4V 


2.0el0 


15.9 


-1.7 


B93 














3.0el0 


20.9 


3.3 


B93 


108767 


5 CrV 


295.5 


46.1 


26.9 


B9.5V 


3.4el0 


-0.5 


7.9 


W96 


22468 


HR 1099 


184.9 


-41.6 


29.0 


G5IV/K1IV 


7.9el7 


21.9 


-2.3 


P97 














1.6el7 


8.2 


-16.0 


P97 














4.0el7 


14.8 


-9.4 


P97 


161868 


7 Oph 


28.0 


15.0 


29.1 


A0V 


3.0el0 


-33.1 


-6.0 


CCW97 














4.2el0 


-29.9 


-2.8 


CCW97 
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Table 3 — Continued 



HD 


Name 


1 


b 


d 


Spec. 


N(X) a 


Velocity 


dVl(96) 


Ref. b 






deer 


deer 


DC 






km s ^ 


km s — ^ 




4128 


(3 Cet 


111.3 


-80.7 


29.4 


KOIII 


5.9el6 


1. 


-5.3 


P97 














2.23el8 


8. 


1.7 


P97 


120418 


9 Peg 


67.4 


-38.7 


29.6 


A2IV 


3.0el0 


-4.2 


4.6 


CLW98 


358 


a And 


111.7 


-32.8 


29.8 


B8IVp 


1.99el0 


13.0 


6.2 


CLW98 


8538 


S Cas 


127.2 


-2.4 


30.5 


A5III-IV 


0.51el0 


12.9 


1.1 


LB92 


120315 


7] UMa 


100.7 


65.3 


30.9 


B3V 


0.9el0 


-3.4 


2.1 


FW02 


209952 c 


a Gru 


350.0 


-52.4 


31.1 


B7IV 


2.0el0 


13.0 


-1.9 


CD95 














1.12el0 


10.2 


0.9 


CD95 


213558 


q Lac 


101.3 


-6.6 


31.4 


A1V 


l.OelO 


3.5 


3.4 


LB92 


88955 


HR 4023 


274.3 


11.9 


31.5 


A2V 


2.38el0 


-1.7 


-4.4 


CLW98 














3.98el0 


10.5 


7.8 


CLW98 


112413 


q2 CVn 


118.3 


78.8 


33.8 


AOspe 


0.51el0 


-1.9 


2.2 


V93 


177756 


A Aql 


30.3 


-5.5 


38.4 


B9Vn 


0.41el0 


-21.9 


3.7 


FW02 














0.41el0 


-26.5 


-0.9 


FW02 














0.96el0 


-30.7 


-5.1 


FW02 


215789 


e Gru 


338.3 


-56.5 


39.7 


A3V 


1.58el0 


-12.2 


-4.3 


CCW97 














1.95el0 


-7.2 


0.7 


CCW97 














2.88el0 


-1.1 


6.8 


CCW97 


218045 


a Peg 


88.3 


-40.4 


42.8 


B9III 


0.8el0 


-4.7 


-3.2 


FW02 














l.lelO 


-0.5 


1.0 


FW02 














0.6el0 


2.5 


4.0 


FW02 


10144 c < e 


a Eri 


290.8 


-58.8 


44.0 


B3Vpe 


0.50el0 


18.9 


16.2 


CD95 














1.26el0 


21.2 


18.5 


CD95 














1.86el0 


11.0 


8.3 


CD95 














2.51el0 


7.6 


4.9 


CD95 


141003 


P Ser 


26.0 


47.9 


46.9 


A2IV 


l.OelO 


-16.7 


5.5 


FW02 














2.9el0 


-23.3 


-1.1 


FW02 














6.6el0 


-20.7 


1.5 


FW02 


212061 


7 Aqr 


62.2 


-45.8 


48.4 


AOV 


7.94e9 


-4.5 


3.9 


CCW97 


135742 


/? Lib 


352.0 


39.2 


49.1 


B8V 


0.8el0 


-33.7 


-10.1 


FW02 














6.4el0 


-23.6 


0.0 


FW02 
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Table 3 — Continued 



HD 


Name 


1 


b 


d 


Spec. 


N(X) a 


Velocity 


dCi(96) 


Ref. b 






deg. 


deg. 


pc 






km s -1 


km s _1 
















6.4el0 


-26.9 


-3.3 


FW02 


106625 


7 CrV 


291.0 


44.5 


50.0 


B8IIIp 


1.5el0 


1.6 


8.5 


W96 














5.4el0 


-2.0 


4.9 


W96 


148857 


A Oph 


17.1 


31.8 


50.9 


AOVpl 


4.7el0 


-24.8 


1.5 


FW02 


160613 


o Ser 


13.3 


9.2 


51.5 


A2Va 


5.6el0 


-29.0 


-0.9 


FW02 


222439 


k And 


109.8 


-16.7 


52.0 


B9IVn 


l.lelO 


0.8 


-4.2 


V93 














6.1el0 


7.6 


2.6 


V93 


186882 


S Cyg 


78.7 


10.2 


52.4 


B9.5IV 


2.8el0 


-18.8 


-6.9 


W96 














2.9el0 


-16.3 


-4.4 


W96 














3.2el0 


-9.6 


2.3 


W96 


62044 


a Gem 


191.2 


23.3 


55.5 


K1III 


5.85el7 


32.0 


8.9 


D97 














7.50el7 


21.4 


-1.7 


D97 


193924 


a Pav 


340.9 


-32.5 


56.2 


B2IV 


1.24el0 


-19.6 


-2.8 


CLW98 














1.05ell 


-18.6 


-1.8 


CLW98 


213998 


r\ Aqr 


66.8 


-47.6 


56.3 


B9IV-Vn 


2.70el0 


-2.1 


4.6 


V93 


156928 


v Ser 


10.6 


13.5 


59.3 


A0/A1V 


4.4el0 


-27.7 


0.4 


FW02 


WD0501+527 


G191-B2B 


156.0 


+7.1 


68.8 


DAw 


1.86el8 


19.3 


-2.0 


Sa99 














3.98el7 


8.6 


-12.7 


Sa99 


111812 


31 Com 


114.9 


89.6 


90.9 


G0II 


7.66el7 


-2.4 


3.2 


D97 


52089 


e CMa 


239.8 


-11.3 


132. 


A1V 


4.0el7 


9.2 


-10.2 


GJ01 f 














6.0el7 


16.2 


-3.2 


GJ01 



a Cohmm Densities less than 10 13 are for Ca + , and values greater than this are based on iV(H°) 
measurements. 

b References are listed in Table 2. 

c Star omitted from restricted sample used for determining Vfl w(93), see text. 

d a Eri is the brightest known Be emission line star, with evidence of variable chromosphere and variable 
stellar radial velocities (Porri & Stalio 1988). 

e Star with infrared excess indicating dust disk (e.g. Habing et al. 2001). 
f Velocities are from GJ01 text and not GJ01 table. 
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Table 4. Regional Properties of Flow 





Cloud a 


Vector 


dVJ Range 


Location 15 


Stars Showing 




Name 




(km s _1 ) 


d ■ l ll b n 


Cloud Component 




1 


LIC 


WHe) 


0.0±1.5 


pc, Z n =150°^250° 


All stars in interval 


T 


Blue Cloud 


Vfiow(96) 


-9.8±0.3 


<3 pc, 233 o ±7°,10 o ±2° 


a CMa, e CMa 


3 


Aql/Oph 


Vfiow(96) 


-6.1±1.0 


<5 pc, 38 o ±10°,9°±14 o 


a Aql, a Oph, ( Aql, 7 Oph, 
A Aql (5 Cyg) 


4 


Peg/Aqr 


4>w(96) 


4.2±0.5 


<30 pc, 75 ±13°,-44 ±5° 


9 Peg, ct Peg, 7 Aqr, 77 Aqr 


5 


North Pole 


W96) 


1.7±0.6 


<22 pc, 6 n >+47°, 


7 UMa, a CrB, 5 UMa, r] UMa, 
a 2 CVn, /3 Ser 


6 


South Pole: 


V flow (96) 


2.4±1.0 


<3.6 pc, 6 n <-48° 


e Ind, a Hyi, r 3 Eri, /3 Cet 


7 


Filament: 


W96) 


6.7±1.8 


<22 pc, 284 o ±12°,-4°±50 o 


5 CrV, 7 CrV, HR4023, 5 Vel, 
a Hyi, e Gru 



a Clouds marked with a ":" are poorly identified and uncertain. 
h d min is the minimum cloud distance. 

c Cloud 2 is the "Blue Cloud" (Lallement et al. 1994; Gry & Jenkins 2001). 



-27- 



Table 5. LISM Heliocentric Velocity Vectors 



Vector 


V 


1 


b 


Basis 


Rcf. 




km s" 1 


o 


o 






Crutcher (C) 


-28±2 


25 


+10 


Optical Ti + , d<100 pc stars 


Crutcher (1982) 


Frisch & York (FY) 


-27 


34 


+15 


Optical Ca + , d<100 pc stars 


Frisch & York (1986) 


LIC 


26+1 


6+3 


16+3 


Ca+ 


LB92, Bertin et al. (1993) 


G 


29 


4.5 


20.5 


Ca+ 


Bertin et al. (1993) 


Frisch (F) 


-26.8±1.3 


6.2 


11.7 


Ca + , d<50 pc stars 


Frisch (1995) 


Bulk Flow (BF) 


-28.1±4.6 


12.4 


11.6 


UV, Optical 


This work 



J I I I I I I I I I I I I L 

20 40 60 

Distance (pc) 



Fig. 1. — (a) Dispersion (dVi(96)) of component velocities in the rest frame of the CLIC as 
defined by Vq ow (96) (see Table 1). (b) Plot of dV[(He) versus star distance, for components 
within 1.5 km s^ 1 of V fl0 w(He) (i.e. |dVS(He)| <1.5 km s" 1 ). 
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Fig. 2. — (a) The galactic coordinates of the background stars used to derive Vf[ ow (96) (Table 
1)1 (b) Same stars as in (a), but coded for dV;(96): circle - |dVi(96)| < 5 km s -1 ; box - 
dV;(96) > 5 km s -1 ; cross - dVj(96) < -5 km s^ 1 . Stars with multiple components show 
more than one symbol, (c) The stars for which |dl / ;(96)| <5 km s _1 are coded for the sign 
of dVi(96): box - dM(96)>0 km s" 1 ; cross - d\X(96)<0 km s _1 . Stars with components in 
each range show both symbols, (d) The galactic coordinates of the 30 stars with components 
showing |dVi(He)| <1.5 km s _1 . Boxes surround the positions of stars within 15 pc. 
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Fig. 3. — (a) dVi(96) plotted against galactic longitude for the restricted sample of stars, 
(b) dVi(LSR) plotted against galactic longitude for same stars. The sinusoidal appearance 
shows that the ISM within 30 pc of the Sun is not near the rest velocity of the LSR (e.g. 
also see Frisch (1995)). (c) dV^(He) plotted against galactic longitude for same stars. 
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Fig. 4. — Histogram of devalues, calculated used Vfl ow (96). A best-fit exponential compo- 
nent distribution (dotted line) is overplotted for r]—5. 1 (Equation 5). The best-fit Gaussian 
distribution (solid line) for these same components, with 6=6.2 km s -1 , is also shown (Equa- 
tion 4). 



Fig. 5.— The V Row (96) upstream direction in the LSR (Z n =1.8°, 6 n =-3.2°) is shown plotted 
on a map of 21 cm emission integrated over the velocity range -450 to 400 km s" 1 , with the 
H° Loop I faintly visible (Hartmann & Burton 1997). The Loop I position in 408 MHz radio 
continuum emission is overplotted as two lines (from Berkhuijsen 1971; Haslam et al. 1982). 
The LSR upstream direction of Vq ow (96) (dot) is ~20° from the tangential direction of Loop 
I, consistent with a CLIC origin associated with either an expanding Loop I superbubble, or 
an outflow from the Sco-Cen Association (e.g. Frisch 1995). 
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Fig. 6. — The Clouds listed in Table 4 are shown plotted against positions of stars in the 
restricted sample. Positions of stars showing Cloud 1 components are circled, and other 
clouds are identified by cloud number. The LSR upstream direction of Vg ow (96) is identified 
by the circled cross. 
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Fig. 7. — Heliocentric upstream directions for various bulk velocity vectors found for nearby 
ISM are plotted. Labels and data for individual vectors are presented in Table 5. 
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